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ABSTRACT 

Massive star formation exhibits an extremely rich chemistry. However, not much evolutionary details 
are known yet, especially at high spatial resolution. Therefore, we synthesize previously published 
Submillimeter Array high-spatial-resolution spectral line observations toward four regions of high- 
mass star formation that are in various evolutionary stages with a range of luminosities. Estimating 
column densities and comparing the spatially resolved molecular emission allows us to characterize 
the chemical evolution in more detail. Furthermore, we model the chemical evolution of massive 
warm molecular cores to be directly compared with the data. The four regions reveal many different 
characteristics. While some of them, e.g., the detection rate of CH3OH, can be explained by variations 
of the average gas temperatures, other features are attributed to chemical effects. For example, C'^'^S 
is observed mainly at the core-edges and not toward their centers because of temperature-selective 
desorption and successive gas-phase chemistry reactions. Most nitrogen-bearing molecules are only 
found toward the hot molecular cores and not the earlier evolutionary stages, indicating that the 
formation and excitation of such complex nitrogen-bearing molecules needs significant heating and 
time to be fully developed. Furthermore, we discuss the observational difhculties to study massive 
accretion disks in the young deeply embedded phase of massive star formation. The general potential 
and limitations of such kind of dataset are discussed, and future directions are outlined. The analysis 
and modeling of this source sample reveals many interesting features toward a chemical evolutionary 
sequence. However, it is only an early step, and many observational and theoretical challenges in that 
field lie ahead. 

Subject headings: stars: formation - stars: early-type - stars: individual (Orion-KL, G29.96, 
IRAS 23151-^5912, IRAS 05358-1-3543) - ISM: molecules - ISM: lines and bands 
- ISM: evolution 



1. INTRODUCTION 

One of the main interests in high-mass star formation 
research today is a thorough charac terization of the ex- 
pectcd massive accretion disks (e.g., Yorke fc Son nhalteil 
[2002: Krumholz et al.l [2007,) . Observations of this type 
require high spatial resolution, and disks have been 
suggested in numerous systems, although not always 
with the same trac ers (for a recent compilation see 
ICesaroni et al.l l2007f ). A major complication is the 
large degree of chemical diversity in these regions. This 
is seen on large scales in terms of the large num- 
ber of molecular detections (e.g., ISchilke et al.l 119971 : 
Ivan Dishoeck fc Blakd 119981 ) but in particular on the 
small scale (of the order 10"* AU and even smaller) 
where significant che mical differentiation is found (e.g., 
iBeuther et al.|[2005al) . At present, we are only beginning 
to decipher the chemical structure of these objects at 
high spatial resolution, whereas our understanding of 
the physical properties better allo ws us to place regions 
into an evolutionary con text (e.g.. IBeuther et ani2007£3 : 
iZinnecker fc Yorki 120071 ). It is not yet clear how our 
emerging picture of the physical evolution is related to 
the observed chemical diversity and evolution. 

Various theory groups work on the chemica l evolu- 
tion during massive star formation (e.g., Casel li et al.l 
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IT993'; 'Millar ct al."1997'; 'CharnleV'1997"; 'Viti ct al."200i 
i Nom ura & Millar 2004; Wakelam e t al. 2005; Dotv et all 
I2OO2I [20061 ). and the results are promising. However, 
the observational database to test these models against 
is still relatively poor. Some single-dish low-spatial- 
resolution line surveys toward several sources do ex- 
ist, but they are all conducted with different spatial 
resolution and covering different frequency bands (e.g., 
Blake e t al.' '19 871 iMacDonald et all 119961: ISchilke et all 



19971: Hatche h et al.l [19981: iMcCutcheon et al.l 1 20001: 
^n der T ak et all l2000l l2003t iJohnstone et all [2003 : 
Bisschop et al.l [20071 ). Furthermore, the chemical struc- 
ture in massive star-forming regions is far from uni- 
form, and at high resolution one observes spatial 
variations between many species, prominent examples 
are Orion- KL, W3OH/H2O or Ccpheus A (see , e.g., 
Wright et al. 1996: W vrowski ct al. 1999; Beuth er et all 
2005aHBrogan et al.l^2007^ . 



Single-dish studies targeted larger source samples at 
low spatial resolution and described the averaged chem- 
ical p roperties of the target regions fe.g.. iHatchell et al] 
Il998t Bisschop et al.l l2007t ). However, no consistent 
chemical investigation of a sample of massive star- 
forming regions exists at high spatial resolution. To ob- 
tain an observational census of the chemical evolution at 
high spatial resolution and to build up a database for 
chemical evolutionary models of massive star formation, 
it is important to establish a rather uniformly selected 
sample of massive star-forming regions in various evo- 
lutionary stages. Furthermore, this sample should be 
observed in the same spectral setup at high spatial reso- 
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lution. While the former is necessary for a rehable com- 
parison, the latter is crucial to disentangle the chemical 
spatial variations in the complex massive star-forming re- 
gions. Because submm interferometric imaging is a time 
consuming task, it is impossible to observe a large sam- 
ple in a short time. Hence, it is useful to employ synergy 
effects and observe various sources over a few years in 
the same spectral lines. 

We have undertaken such a chemical survey of mas- 
sive molecular cores containing high-mass protostars in 
different evol utionary s tages using the Submillimeter Ar- 
ray (SMA*, iHo et all 12004) since 2003 in exactly the 
same spectral setup. The four massive star-forming re- 
gions span a range of evolutionary stages and luminosi- 
ties: (1) t he prototypical hot niolecul ar core (HMC) 
Orion-KL (iBeuther et a ll l2004bL [20051. (2'! an HMC 
at larger distance G29.96 (|Beuther et al.l l2007cD . and 
two regions in a presumably earlier evolutionary phase, 
namely (3) the younger but similar luminous High- 
Mass Prot ostcUar Object (HMPO) IRAS 23151+5912 
(jBeuther et al. 200 ?d) and (4) the less luminous HMPO 
IRAS 05358-^3543 (jLeurini et all I2007D . Although the 
latter two regions also have central temperatures >100K, 
qualifying them as "hot" , their molecular line emission is 
considerably weaker than from regions which are usually 
termed HMCs. Therefore, we refer to them from now 
o n as early-HMPOs (se e also the evolutionary sequence 
in IBeuther et al.ll2007a| ). Table [1] lists the main physical 
parameters of the selected target regions. 

The SMA offers high spatial resolution (of the order 1") 
and a large enough instantaneous bandwidth of 4 GHz 
to sample numerous molecular transitions simultaneously 
(e.g., ^^SiO and its rarer isotopologue '^"SiO, a large se- 
ries of CH3OH lines in the Vt = 0,1,2 states, CH3CN, 
HCOOCH3, SO, SO2, and many more lines in the given 
setup, see ^J2]). Each of the these observations have been 
published separately where we provide d etailed discus- 
sions of the particularities of each sourc e (jBeuther et al.l 
l2004bl[2005all2007dHLeurini et al.l[200l V These objects 
span an evolutionary range where the molecular gas and 
ice coated grains in close proximity to the forming star 
are subject to increasing degrees of heating. In this fash- 
ion, volatiles will be released from ices near the most 
evolved (luminous) sources altering the surrounding gas- 
phase chemical equilibrium and molecular emission. In 
this paper, we synthesize these data and effectively re- 
observe these systems at identical physical resolution in 
order to identify coherent trends. Our main goals are 
to search for (a) trends in the chemistry as a function 
of evolutionary state and (b) to explore, in a unbiased 
manner, the capability of molecular emission to trace co- 
herent velocity structures that have, in the past, been 
attribute to Keplerian disks. We will demonstrate that 
the ability of various tracers to probe the innermost re- 
gion, where a disk should reside, changes as a function 
of evolution. 

2. DATA 

The four sources were observed with the SMA between 
2003 and 2005 in several array configurations achieving 

The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



(sub)arcsecond spatial resolution. Fo r detailed obs er- 
vati onal description, see IBeuther et al.l (|2005al l2007cl ldl) 
and lLeurini et al.l (|2007D . The main point of interest to 
be mentioned here is that all four regions were observed 
in exactly the same spectral setup. The receivers op- 
erated in a double-sideband mode with an IF band of 
4-6 GHz so that the upper and lower sideband were sep- 
arated by 10 GHz. The central frequencies of the upper 
and lower sideband were 348.2 and 338.2 GHz, respec- 
tively. The correlator had a bandwidth of 2 GHz and 
the channel spacing was 0.8125 MHz, resulting in a nom- 
inal spectral resolution of ~0.7kms~^. However, for the 
analysis presented below, we smoothed the data-cubes to 
2kms~^. The spatial resolution of the several datasets 
is given in Table [2] Line identifications were done in an 
iterative way: we first com pared our data with t he single- 
dish line survey of Orion bv lSchilke et al.l (|1997f ) and then 
refined the analysis via the molecular spectroscopy cata- 
logs of JPL and t he Cologne database of molecular spec- 
troscopy CDMS (jPovnter fc Pickettlfl985l : iMueller et all 
,2002) . 

3. RESULTS AND DISCUSSION 

3.1. Submm continuum emission 

To set the four regions spatially into context. Figure 
[1] presents the four submm continuum images obtained 
with the SMA. As one expects in a clustered mode of 
massive star formation, all regions exhibit multiple struc- 
tures with the number of sub-sources > 2. As discussed 
in the corresponding papers, while most of the submm 
continuum peaks likely correspond to a embedded pro- 
tostellar sources, this is not necessarily the case for all of 
them. Some of the submm continuum peaks could be ex- 
ternally heated gas cl umps (e.g., the Orion hot core peak, 
IBeuther et aIll2004bD or may be produced by shock inter- 
actions with inolecula r outflows (e.g., IRAS 05358-1-3543, 
IBeuther et al.ll20"07b| ). In the following spectral line im- 
ages, we will always show the corresponding submm con- 
tinuum map in grey-scale as a reference frame. 

3.2. Spectral characteristics as a function of evolution 

Because the four regions are at different distances (Ta- 
ble [1]), to compare the overall spectra we smoothed 
all datasets to the same linear spatial resolution of 
~5700AU. Figure 121 presents the final spectra extracted 
at this common resolution toward the peak positions of 
all four regions. Furthermore, we present images at the 
original spatial resolution of the four molecular species or 
vibrationally-torsionally excited lines that are detected 
toward all four target regions (Figs. 3 to 6). 

The spectral characteristics between the HMCs and 
the early-HMPOs vary considerably. Table |4| lists all de- 
tected lines in the four regions with their upper energy 
levels Eu/k and peak intensities S'poak at a spatial reso- 
lution of -5700 AU. 

3.2.1. Excitation and optical depth effects? 

Since the line detections and intensities are not only 
affected by the chemistry but also by excitation effects, 
we have to estimate quantitatively how much the latter 
can influence our data. In local thermodynamic equilib- 
rium, the line intensities are to first order depending on 

the Boltzniann factor e~i^ and the partition function Q: 
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TABLE 1 

Source parameters 





Orion-KL^' 


G29.96'' 


23151'' 


05358'' 


L- [Lq] 


10^ 


9 X 10" 


10^ 


103.8 


d[pc] 


450 


6000 


5700 


1800 


[Mq] 




2500"= 


600 


300 


tL [k] 


300 


340 


150 


220 


A^pcak(H2)Mcm-2] 


9 X 1024 


6 X 1024 


2 X 10^"' 


2 X lO^-* 


Type 


HMC 


HMC 


early-HMPO 


early-HMPO 



" Luminosities and masses are derived from single-dish data. Since most regions split up into multiple sources, individual values for sub 

members are lower. 

' The SMA da ta are first published in Eeuther et al. 2005a, 2007c d and Leurinietal] ll2007h 
[Menten & Reldl | |T995|) : lOlmi et al.l | |2003) : ISridharaiTet al. (2002); Bcutlicr et al. (20 02b,) . 



Other parameters are taken from 



TABLE 2 
Spatial resolution 



The integrated masses should be accurate within a factor 5 (^cuthcr ct al. 2002b). 
This value was calculated from the 870 lira flux of lSchilke eFal.. (1997 ) following .Hildebrandl lll983l) assuming an average temperature of 
50K. 

*= This value was calculated from the 850 /.tm flux of lThompson et al.l ll200l) following iHildebrandl 1198^ as in comment b. 

f Peak rotational temperatures derived from CH3OH. 

s H2 column densities toward the peak positions derived from the submm dust continuum observatios l IBeuther et al.ir2004bl . l2007d ldllH) . 

As shown in Table IH our spectral setup barely con- 
tains few lines from rarer isotopologues. We therefore 
cannot readily determine the optical depth of the molec- 
ular lines. While it is likely that, for example, the ground 
state CH3OH lines have significant optical depths, rarer 
species and vibrationally excited lines should be more op- 
tically thin. We checked this for a few respective species 
(e.g., C'^'^S or SO2) via running large -velocity gradient 
models (LVG, Ivan der Tak et al1l2007t l with typical pa- 
rameters for these kind of regions, confirming the over- 
all validity of optically thin emission for most lines (see 
^3.2.21 fc i j3.2.4p . However, without additional data, we 
cannot address this issue in more detail. 





Orion-KL 


G29.96 


23151 


05358 


Cont. ["] 


0.8 X 0.7 


0.4 X 0.3 


0.6 X 0.5 


1.1 X 0.6 


Av. Cont. [AU] 


340 


2100 


3100 


1500 


Line ["] 


1.4 X 1.1 


0.6 X 0.5 


1.1 X 0.8 


1.4 X 0.8 


Av. Line [AU] 


560 


3300 


5400 


2000 



I{T)dv (X 



Q{t) 



where is the upper level energy state. For poly- 
atomic molecules in the high-temperature limit, one can 
approximate Q{T) cx VT^ (e.s.. iBlake et al.l ll98/). To 
get a feeling how much temperature changes affect lines 
with different one can form the ratio of / I{T)dv at 
two different temperatures: 



JI{T)dv 

Equating this ratio now for a few respective upper level 
energy states and gas temperatures of T = 100 & 50K 
(Table [3]), we find that the induced intensity changes 
mostly barely exceed a factor 2. Only for very highly 
excited lines like CH30H(7ij — 6ifi)vt = 1 at relatively 
low temperatures (T = 50 K) do excitation effects be- 
come significant. However, at such low temperatures, 
these highly excited lines emit well below our detection 
limits, hence this case is not important for the present 
comparison. Therefore, excitation plays only a minor 
role in producing the molecular line differences discussed 
below, and other effects like the chemistry turn out to be 
far more important. 



TABLE 3 

Excitation effects 



Line 


-Ell 
k 

(K) 


J I(2T)dv 
J I{T)dv 

@ T=100K 


J I{2T)dv 
J I{T)dv 

@ r=50K 


C=*''S(7-6) 


65 


0.49 


0.68 


302(184,1-183,1) 


197 


0.95 


2.5 


CH30H(7i,7 -6i,6)t)t = 1 


356 


2.1 


12.4 



3.2.2. Column densities 

Estimating reliable molecular column densities and/or 
abundances is a relatively difficult task for interferomet- 
ric datasets like those presented here. The data are suf- 
fering from missing short spacings and filter out large 
fractions of the gas and dust emission. Because of the dif- 
ferent nature of the sources and their varying distances, 
the spatial filtering affects each dataset in a different way. 
Furthermore, because of spatial variations between the 
molecular gas distributions and the dust emission rep- 
resenting the H2 column densities, the spatial filtering 
affects the dust continuum and the spectral line emis- 
sion differently. On top of this, Figures [3] to [6] show that 
is several cases the molecular line and dust continuum 
emission are even spatially offset, preventing the estima- 
tion of reliable abundances. 

While these problems make direct abundance estimates 
relative to H2 impossible, nevertheless, we are at least 
able to estimate approximate molecular column densities 
for the sources. Since we are dealing with high-density 
regions, we can derive the column densities from the spec- 
tra shown in Figure [2] assuming local thermodynamic 
equilibrium (LTE) and optically thin emission. We mod- 
eled the molecular emission of each species separately 
using the XCLASS superset to the CLASS software de- 
veloped by Peter Schilke (priv. comm.). This software 
package uses the l ine c atalogs from JPL and CDMS 
(|Povnter fc Pickettlll985l : iMiiller et al.ll200lD . The main 
free parameters for the molecular spectra are tempera- 
ture, source size and column density. We used the tem- 
peratures given in Table [1] (except of Orion-KL where we 
used 200 K because of the large smoothing to 5700 AU) 



4 



Beuther et al. 



Orion K 



o 

CN 
CN 

o 

LD 

I 




0.4 



0,2 



CN 
o 

CN 



G29.96 



o 



- 



lO 

CN 
CN 

o 

I 5 35 1 4.6 

IRAS23151 +5912 



-0.2 



CN 
CJ) 

o 

CN 




6000AU 



] 



0.1 5 







- 0.1 



0.05 



- 



1 4'4 



18^46"'03!8 03'7 
RAS05358+3543 



00 

CO 
CN 
o 

CJ5 

LO 



CD 

00 
CN 
o 

CD 
LO 



o 

6000AU 




0.05 







LO 
LD 

LO 
■+ 
o 

LO 



O 

LO 
LO 
o 

LO 



6000AU 





5^39"'l 3'0 



o 



1 2'5 



0.1 5 



0. 



0.05 







Fig. 1. — 862 /im continuum images toward tiie four massive star-forming regions. Positive and negative features are shown as solid 
and dashed contours, respectively. The contouring always starts at the 3cr levels and continues in 2<j steps (with la values of 35, 21, 5, 
7mjybeam~^, respectively). These figures are adaptions from the papers by Beut her etaTI ([gOOSa, 2007c, d b). Scale-bars are presented 
in each panel, and the wedges show the flux density scales in Jy. The axis are labeled in R.A. (J2000) and Dec. {J2000). The synthesized 
beams sizes are listed in Table [2] 



with approximate source sizes estiraated from the dust 
continuum and spectral hue maps. Then we produced 
model spectra with the column density as the remain- 
ing free parameter. Considering the missing flux and 
the uncertainties for temperatures and source sizes, the 
derived column densities should be taken with caution 
and only be considered as order-of-magnitude estimates. 
Table [5] presents the results for all sources and detected 
molecules. 

3.2.3. General differences 

An obvious difference between the four regions is 
the large line forest observed toward the two HMCs 
Orion-KL and G29.96 and the progressively less de- 
tected molecular lines toward IRAS 23151+5912 and 
IRAS 05358+3543. Especially prominent is the difference 
in vibrationally-torsionally excited CII3OII lines: we de- 
tect many transitions in the HMCs and only a single 
one in the two early-HMPOs. Since the vibrationally- 
torsionally excited CII3OH lines have higher excitation 



levels Eu/k, this can be relatively easily explained by 
on average lower temperatures of the molecular gas in 
the early-HMPOs. Assuming an evolutionary sequence, 
we anticipate that the two early-HMPOs will eventually 
develop similar line forests like the two HMCs. 

Analyzing the spatial distribution of CH3OH we find 
that it is associated with several physical entities (Figs. [3] 
& [4]). While it shows strong emission toward most 
submm continuum peaks, it exhibits additional inter- 
esting features. For example, the double-peaked struc- 
ture in G29.96 (Fig [3]) may be caused by high optical 
depth of the molecular emission, whereas the lower opti- 
cal depth vibrationally-torsionally excited lines do peak 
toward the central dust and gas core (Fig 21 smoothing 
the submm continuum map to the lower spatial resolu- 
tion of the line data, the four submm sources merge into 
one central peak). In contrast, toward Orion-KL the 
strongest CH3OH features in the ground state and the 
vibrationally-torsionally excited states are toward the 
south-western region called the compact ridge. This is 
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Fig. 2. — SMA spectra extracted from the final data-cubes in the image domain toward four massive star-forming regions (Orion-KL 
top-row, G29.96 second row, IRAS 23151+5912 third row and IRAS 05358-1-3543 fourth row). For a better comparison, the data-cubes 
were all smoothed to the same spatial resolution of ~5700AU from the current IRAS 23151+3543 dataset. The circles in Fig. |3] outline the 
corresponding spatial regions. The spectral resolution in all spectra is 2km/s. 
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the interface between a molecular outflow and the am- 
bient gas. Our data confirm previo us work which sug - 
gests an abundance enrichment fe.g. jBlake et al.lll987D . 
For instance in qui escent gas in Orion t he ratio of 
CHgOH/C^^S <20 (jBergin fc Langeii fTool . while our 
data find a ratio of 100 (see Table [5]). This is believed 
to be caused by outf low shock processe s in the dense 
surrounding gas (e.g.. IWright et allll996[ ). Furthermore, 
as will be discussed in §3.3( there exist observational in- 
dications in two of the sources (IRAS 23151-1-5912 and 
IRAS 05358-f 3543) that the vibrationally-torsionally ex- 
cited CH3OII may be a suitable tracer of inner rotating 
disk or tori structures. 

Toward G29.96, we detect most lines previously also 
observed toward Orion-KL, a few exceptions are the 
■^°SiO line, some of the vibrationally-torsionally state 
vt — 2 CH3OII lines, some N-bearing molecular lines 
from larger molecules fike CH3CH2CN or CH3CHCN, 
as weU as a few ^^S02 and HCOOCH3 lines. While 
for some of the weaker lines this difference may partly 
be attributed to the larger distance of G29.96, for other 
lines such an argument is unlikely to hold. For example, 
the CII3CII2CN line at 348.55 GHz is stronger than the 
neighboring II2CS line in Orion-KL whereas it remains 
undetected in G29.96 compared to the strong II2CS line 
there. This is also refiected in the different abundance 
ratio of CH3CII2CN/II2CS which is more than an order 
of magnitude larger in Orion-KL compared with G29.96 
(see TableEJ. Therefore, these differences are likely trac- 
ing true chemical variations between sources. In contrast, 
the only lines observed toward G29.96 but not detected 
toward Orion-KL are a few CII3OCH3 lines. 

The main SiO isotopologue ^^810(8-7) is detected in 
aU sources but IRAS 053584-3543. This is relatively 
surprising because SiO(2-l) is strong in this region 
(jBeuthcr ct al. 2GG2ai), and the upper level energy of 
the J = 8 — 7 of ~75K does not seem that extraor- 
dinarily high to produce a non-detection. For exam- 
ple, the detected CH30II(7ij — 6ifi)vt — 1 transition 
has an upper level energy of 356 K. This implies that 
IRAS 05358-1-3543 does have warm molecular gas close 
to the central sources. However, the outfiow-components 
traced by SiO are at on average lo wer temperatures 
(prob ably of the order 30 K, e.g., ICabrit fc Bertout) 
fl990h which may be the cause of the non-detection in 
IRAS 05358-f 3543. Furthermore, the critical density of 
the SiO (8-7) line is about two orders of magnitude higher 
than that of the (2-1) transition. Hence, the density 
structure of the core may cause the (8-7) non-detection 
in IRAS 05358-^3543 as well. 

While the rarer ■^*'SiO(8-7) isotopologue is detected to- 
ward Orion-KL with nearly co mparable stren gth as the 
main isotopologue (Fig. [2l and .Beuther et al.ir2005a ). we 
do not detect it at all in any of the other sources. 

A little bit surprising, the H2CS line at 338.081GHz 
is detected toward Orion-KL, G29.96 as well as the 
lowest luminosity source IRAS 05358-f 3543, however, it 
remains undetected toward the more luminous HMPO 
IRAS 23151-f 5912. We are currently lacking a good ex- 
planation for this phenomenon because H2CS is pre- 
dicted by most chemistry networks as a parent molecule 
to be found early in t he evolutionary sequence (e.g., 
iNomura fc Millarl l2004f ). The sulphur and nitrogen 
chemistries are also peculiar in this sample, and we out- 



line some examples below. 

3.2.4. Sulphur chemistry 

The rare Carbon-sulphur isotopologue C'^^S is detected 
toward all four regions (Fig. [2|). However, as shown in 
Fig.[5]C^''S does not peak toward the main submm con- 
tinuum peaks but is offset at the edge of the core. In the 
cases of G29.96 and IRAS 23151-f 5912 the C^^S morphol- 
ogy appears to wrap around the main submm continuum 
peaks. Toward IRAS 05358-f 3543 C^^S is also weak to- 
ward the strongest submm peak (at the eastern edge of 
the image) but shows the strongest C^'^S emission fea- 
tures offset from a secondary submm continuum source 
(in the middle of the image). Toward Orion-KL, weak 
C'^^'S emission is detected in the vicinity of the hot core 
peak, whereas we find strong C'^^S emission peaks offset 
from the dust continuum emission. 

To check whether our optical-thin assumption from 
i 33.2.2l is valid also for C^'^S we ran LVG rad iative transfer 
models (RADEX. I^an der Tak et al.ll2007f ). We started 
with the H2 column densities from Table [T] and assumed 
a typical CS/H2 abu ndance of 10~^ with a terrestrial 
CS/C^^S ratio of 23 (| Wannieii [l98l . Above the criti- 
cal density of 2 x lO^cm"^, with the given broad C'^^S 
spectral FWHM (between 5 and 12kms^"'^ for the four 
sources), the 0^4(7-6) emission is indeed optically thin. 
Hence optical depth effects are not causing these large 
offsets. Since furthermore the line intensities depend not 
just on the excitation but more strongly on the gas col- 
umn densities, excitation effects only, as quantified in 
t ^3.2.11 cannot cause the observational offsets as well. 
Therefore, chemical evolution may be more important. 
A likely scenario is based on different desor ption tem- 
peratures of molecules from dust grains (e.g.. IViti et al.l 
i2004) : CS and C^'^S are desorbed from grains at tem- 
peratures of a few 10 K, and at such temperatures, these 
molecules are expected to be well correlated with the dust 
continuum emission. Warming up further, at 100 K H2O 
desorbs and then dissociates to OH. The OH quickly re- 
acts with the sulphur forming SO and SO2 which then 
will be centrally peaked (see §3.2.6p . Toward G29.96, 
IRAS 23151-f 5912 and IRAS 05358-f 3543 we find that 
the SO2 emission is centrally peaked toward the main 
submm continuum peaks confirming the above outlined 
chemical scenario (Fig. [5] & !j3.2.6l ') . To further in- 
vestigate the C'^'*S/S02 differences, we produced column 
density ratio maps between C'^'^S and SO2 for G29.96, 
IRAS 23151-f 5912 and IRAS 05358-f 3543, assuming lo- 
cal thermodynamic equilibrium and optically thin emis- 
sion (Fig. [7|). To account for the spatial differences 
that SO2 is observed toward the submm peak positions 
whereas C'^'^S is seen more to the edges of the cores, for 
the column density calculations we assumed the temper- 
atures given in Table [1] for SO2 whereas we used half 
that temperature for C'^'^S. Although the absolute ra- 
tio values are highly uncertain because of the different 
spatial filtering properties of the two molecules, qualita- 
tively as expected, the column density ratio maps have 

^ For G29.96 the SO2 peaks are actually right between the bet- 
ter resolved submm continuum sources. This is due to the lower 
resolution of the line data because smoothing the continuum to the 
same spatial resolution, t hey peak very close t o the SO2 emission 
peaks (see, e.g., Fig. 2 in lBeuther et al.ll2007dV 
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Fig. 3. — CH3OH contour images (line blend between CH30H{72,5 — 62,4) and CH30H(72,6 — 62,5)) toward the four massive star-forming 
regions. Positive and negative features are shown as sohd and dashed contours. The contouring is done from 15 to 95% (step 10%) of the 
peak emission, and the peak emission values are 7.5, 0.9, 1.0 and 0.6 Jybeam"^ from left to right, respectively. The integration regimes 
for the four sources are [5,15], [90,104], [-60,-52] and [-20,-12] kms~^. The grey-scale with dotted contours shows the submm continuum 
emission from Fig. [l] This figure is an adaption from the papers by IBeut her et al. (2005a, 2007c d); Leurini et al. (2007). The axis are 
labeled in R.A. (J2000) and Dec. (J2000). The spatial resolution is listed in Table[2] The circles represent the regions of diameter 5700 AU 
used for the comparison spectra in Fig. [2] 
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Fig. 4. — Rotationally-torsionally excited CH30H(7i,7 — 6i,6)(ft = 1) images toward the four massive star-forming regions. Positive and 
negative features are shown as solid and dashed contours. The contouring is done from 15 to 95% (step 10%) of the peak emission, and 
the peak emission values are 7.5, 0.9, 1.0 and 0.6 Jybeam"^ from left to right, respectively. The integration regimes for the four sources 
are [3,13], [91,105], [-58,-54] and [-18,-10] km s ~^. The grey- scale with dotted contou rs shows the su bmm continuum emission from Fig. [Tl 
This figure is an adaption from the papers bv lBeuthereta l. (2005a, 2007c, d); Leur ini et al.l ||2007|) . The axis are labeled in R.A. (J2000) 
and Dec. (J2000). The spatial resolution is listed in Table [2] 



the lowest values in the vicinity of the submm contin- 
uum sources and show increased emission at the core 
edges. The case is less clear for Orion-KL which shows 
the strongest SO2 emission toward the south-eastern re- 
gion called compact ridge. Since this compact ridge is 
believed to be caused by the interactio n of a molecular 
outflow with the ambient dense gas fe.g. JLiu et al.|[2"002f ) 
and SO2 is known to be enriched by shock interactions 
with outflows, this shock-outflow interaction may dom- 
inate in Orion-KL compared with the above discussed 
C^'*S/S02 scenario. For more details on the chemical 
evolution see the modeling in ^13.2.61 

3.2.5. Nitrogen chemistry 

It is intriguing that we do not detect any nitrogen- 
bearing molecule toward the two younger HIvIPOs 
IRAS 2315H-5912 and IRAS 05358-f 3543 (FigEland Ta- 
ble [4]). This already indicates that the nitrogen- 
chemistry needs warmer gas to initiate or requires more 
time to proceed. To get an idea about the more 
subtle variations of the nitrogen chemistry, one may 
compare some specific line pairs; For example, the 
HN13C/CH3CH2CN line blend (dominated by mi^^G) 
and the SO2 line between 348.3 and 348.4 GHz are of 
similar strength in the HMC Orion-KL (Fig. The 



same is approximately true for the HIVIC G29.96, al- 
though SO2 is relatively speaking a bit weaker there. The 
more interesting differences arise if one contrasts with 
the younger sources. Toward the 10^ Lq early-HlVIFO 
IRAS 23151-^5912, we only detect the SO2 line and the 
HN^^C/CH3CH2CN hne blend remains a non-detection 
in this source. In the lower luminosity early-HlVIPO 
IRAS 05358-1-3543, both lines are not detected, although 
another SO2 line at 338.3 GHz is detected there. 

Judging from these line ratios, one can infer that 
SO2 is relatively easy to excite early-on in the evolu- 
tion of high-mass star- forming regions. Other sulphur- 
bearing molecules like H2CS or CS are released even ear- 
lier from the grains, bu t SO-type molecules are forme d 
quickly (e.g.. ICharnTevI Il997t [N omura fc ]V[illail[200l . 
In contrast to this, the non-detection of spectral lines 
hke the HN13C/CH3CH2CN line blend in the early- 
HlVIPOs indicates that the formation and excitation of 
such nitrogen-bearing species takes place in an evolu- 
tionary more evolved phase. This may either be due to 
molecule-selective temperature-dependent gas-dust des- 
orption processes or chemical network reactions requir- 
ing higher temperatures. Furthermore, simulations of 
chemical networks show that the complex nitrogen chem- 
istry simply needs more time to be activated (e.g., 
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Fig. 5. — C"^''S(7— 6) images toward the four massive star-forming regions. Positive and negative features are shown as soUd and dashed 
contours. The contouring is done from 15 to 95% (step 10%) of the peak emission, and the peak emission values are 7.5, 0.9, 1.0 and 
0.6 Jy beam~^ from left to right, respectively. The integration regimes for the four sources are [-2,14], [92,104], [-58,-51] and [-17,-13] kms~^ . 
The grey-sca le with dotted cont ours shows the subm m continuum emission from Fig. [T] This figure is an adaption from the papers by 
IBeuther et"a l. (2 005a . 20073 ldh: iLeurini et al.l ||2007|'| . The axis are labeled in R.A. (J2000) and Dec. (J2000). The spatial resolution is 
listed in Table [2] 
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Fig. 6. — 302(144^14 — I8345) images toward the four massive star-forming regions. Positive and negative features are shown as solid 
and dashed contours. The contouring is done from 15 to 95% (step 10%) of the peak emission, and the peak emission values are 7.5, 0.9, 
1.0 and 0.6 Jybeam"^ from left to right, respectively. Only for IRAS 05358-f 1732 the contouring starts at the 35% level because of a worse 
signal-to-noise ratio. The integration regimes for the four sources are [0,20], [94,100], [-60,-50] and [-17,-13] kms~ ^. The grey-scale with 
dotted co ntours sh ows the submm continuum emission from Fig. [T] This figure is an adaption from the papers bv IBeuther et al.l H2005al . 
I2007d ldh: ILeurini et al. (2007i). The axis are labeled in R.A. (J2000) and Dec. (J2000). The spatial resolution is listed in Table|2] 



ICharnlev etld] |200H iNomura fc IVIillad I2004D . Recent 
modeling bv'Garro d et alT i 20071 ) indicates that the grad- 
ual switch-on phase of hot molecular cores is an impor- 
tant evolutionary stage to produce complex molecules. 
In this picture, the HlVICs have switched on their heat- 
ing sources earlier and hence had more time to form these 
nitrogen molecules. 

Comparing just the two HJVICs, we find that the 
CH3CH2CN at 348.55 GHz is strong in Orion-KL but 
not detected in G29.96. Since G29.96 also exhibits less 
vibrationally-torsionally excited CH3OH lines, it is likely 
on average still at lower temperatures than Orion-KL and 
may hence have not formed yet all the complex nitrogen- 
bearing molecules present already in Orion-KL. 

3.2.6. Modeling 

To examine the chemical evolution of warm re- 
gions in greater detail we used the chemical model of 
iBergin fc LangeJ (|1997l ). This model includes gas-phase 
chemistry and the freeze-out to and sublimation from 
grain surfaces. The details of this model are provided in 
that paper with the only addition being the inclusion of 
water ice formation on grain surfaces. The binding ener- 
gies we have adopted are for bare silicate grains, except 
for water ice which is assumed to have a binding energy 
appropriate for hydrogen bonding between frozen water 



molecules (jFraser et al.ll200ll) . To explore the chemistry 
of these hot evaporative regions we have run the model 
for IG^yrs with starting conditions at — 10^ cm^'^ 
and Tgas = T(just=20K. Under these conditions most 
gaseous molecules, excluding H2, will freeze onto the 
grain surface and the ice mantle forms, dominated by 
H2O. This timescale (lO^yrs ) is quite short, but is 
longer than the free-fall time at this density and is cho- 
sen as a representative time that gas might spend at very 
high density. After completion we assume that a massive 
star forms and the gas and dust temperature is raised 
such that the ice mantle evaporates and the gas-phase 
chemistry readjusts. We have made one further adjust- 
ment to this model. Our data suggest that HNC (a rep- 
resentative nitrogen-bearing species) is not detected in 
early-HlVIPOs and that the release of this species (or its 
pre-cursor) from the ice occurs during more evolved and 
warmer stages. Laboratory data on ice desorption sug- 
gest that the process is not as simple as generally used 
in models where a giv en molecule evapora tes at its subli- 
mation temperature (jColhngs et al.l[2003 ). Rather some 
species co-desorb with water ice and the key nitrogen- 
bearing species, NH3, falls into this category (see also 
iViti et al] 120041. We have therefore assumed that the 
ammonia evaporates at the same temperature as water 
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Fig. 7. — Column density ratios between C'^'^S over SO2 toward G29 
Tlie contours present the corresponding submm continuum maps as 
into the central region to better show the ratio variations. 



ice. O ur initial abundances are taken from lAikawa et al.l 
(|1996[ ). except we assume that 50% of the nitrogen is 
frozen in the form of NH3 ice. This assumption is con- 
sistent with two sets of observations. First, detections 
of NH3 in ices towards YS O's find abundan ces of ^ 2- 
7% rel ative to H2O ( BottineUi et al.|[2007 ). Dartois et all 
(|2002l ) find a limit of ~ 5% relative to H2O towards other 
YSO's. Our assumed abundance of NH3 ice is 5 xlO~^ 
relative to H2O, assuming an ice abundance of 10~^ (as 
appropriate for cold sources) this provides an ammonia 
ice abundance of ~5%, which is consistent with ice obser- 
vations. Second, high resolution ammonia observations 
often find high abundances of NH3 in the gas phase to- 
wards liot_coreSj_often_as high as 10~^, relative to H2 
(e.g., ICesaroni et"all I1994D . Pure gas phase chemistry 
will have some difficulty making this high abundance in 
the cold phase and we therefore assume it is made on 
grains during cold phases. 

Figure [8] presents our chemical model results at the 
point of star "turn on" where the gas and dust be- 
comes warmer and ices evaporate. Two different models 
were explored. The first labeled as Tdust=70K ("warm 
model" ) is a model which is insufficient to evaporate wa- 
ter ice, while the second, "hot model" withTdust — ^50K, 
evaporates the entire ice mantle. Much of the chemical 
variations amongst the early-HMPO and HMC phases 
is found for CS, SO, SO2, and HNC (note this line is 
blended with CH3CH2CN) and we focus on these species 
in our plots (along with H2O). It is also important to reit- 
erate that due to differences in spatial sampling between 
the different sources as well as between line and con- 
tinuum emission, we cannot derive accurate abundances 
from these data, but rather can attempt to use the mod- 
els to explain trends. For the warm model the main result 
is that the imbalance created in the chemistry by the re- 
lease of most species, with ammonia and water remaining 
as ices on the grains, leads to enhanced production of CS. 
Essentially the removal of oxygen from the gas (exclud- 
ing the O found in CO) allows for rapid CS production 
from the atomic Sulfur frozen on grains during the cold 
phase. Thus, for early-HMPOs, which might not have a 
large fraction of gas above the water sublimation point 
(T ~ 110 K), the ratios of CS to other species are quite 
large. In the hot model, when the temperature can evap- 
orate both II2O and NII3 ice, ratios between the same 
molecules are orders of magnitude lower (Fig.[H]). There 
is a large jump in the water vapor abundance (and NH3, 
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,96, IRAS 23151-1-5912 and IRAS 05358+3543 are shown in grey-scale, 
in shown the previous figures. For IRAS 05358+3543 we zoom only 

which is not shown) between the warm and the hot model 
driving the chemistry into new directions. CS remains in 
the gas, but not with as high abundance as in the warm 
phase and is gradually eroded into SO and ultimately 
S02- Hence, SO2 appears to be a better tracer of more 
evolved stages. In this sense, even the early-HMPOs can 
be considered as relatively evolved, and it will be impor- 
tant to extend similar studies to even earlier evolutionary 
stages. HNC also appears as a brief intermediate product 
of the nitrogen chemistry. 

The above picture is in qualitative agreement with our 
observations: that CS should be a good tracer in early 
evolutionary states even prior to our observed sample, 
but that it is less well suited for more evolved regions 
like those studied here. Other sulfur-bearing molecules, 
in particular SO2, appear to better trace the warmest gas 
near the forming star. HNC, and perhaps other nitrogen- 
bearing species, are better tracers when the gas is warm 
enough to evaporate a significant amount of the ice man- 
tle. 

3.3. Searching for disk signatures 

While the chemical evolution of massive star-forming 
regions is interesting in itself, one also wants to use 
the different characteristics of molecular lines as tools 
to trace various physical processes. While molecules 
like CO or SiO have reg ularly been used to investigate 
molecular outflows fe.g.. lArce et al.ll2007f) . the problem 
to identify the right molecular tracer to investigate accre- 
tion disks in massive star formation is much more severe 
(e.g.. ICesaroni et al]l2007| ). Major observational obsta- 
cles arise from the fact that disk-tracing molecular lines 
are usually often not unambiguously found only in the ac- 
cretion disk, but that other processes can produce such 
line emission as well. For example, molecular lines from 
CN and HCN are high-density tracers and were believed 
to be good candidates to inv estigate disks in embed - 
ded protostellar sources (e.g., lAikawa fc HCTbsil IMilh . 
However, observations revealed that both spectral lines 
are strongly affected by the associated molecular outflow 
and hence difficult t o use for massive disk st udies (CN 
iBeuther et"aII[2?)04aL HCN iZhang et al J [200I . 

As presented in ICesaroni et al.l (|2007 ). various differ- 
ent molecules have in the past been used to investi- 
gate disks/rotational signatures in massive star forma- 
tion (e.g., CH3CN, C^^S, NH3, HCOOCH3, C^^O, Hi^O, 
see also Fig. [9]). The data presented here add three other 
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Fig. 8. — The left and right columns show results from our modeling of the "cold and warm models" , respectively. The two top-panels show 
the abundance ratios between several species important for our observations versus the time. The bottom panels present the abundances 
relative to H2. Time in this figure refers to time after the onset of massive star formation and hence we do not show any evolution during 
the cold starless phase. 

lution is important for these molecules, temperature ef- 
fects have to be taken into account as well. For ex- 
ample, the torsionally excited CH3OH line traces rotat- 
ing motions in IRAS 23151-H5912 arid IRA S 05358-^3543 
(|Beuther et al.ll2007dt Heurini et al.ll2007l Fig. l]) but it 
is weak and difficult to detect in colder and younger 
sources. Therefore, some initial heating is required to 
employ highly excited lines for kinematic studies. In 
addition to these evolutionary effects, o ptical depth is 



potential disk trace rs (RN^^C and HC3N for G29.96, 
iBeuther et al.l[20073 . and torsionally excit ed CH3QH for 
IRAS 23151-H5912 and IR AS 05358+3543 IBeuther et all 
I2007dt iLeurini et all [20071 Fig. E]). An important point 
to note is that in most sources only one or the other spec- 
tral line exclusively allows a study of rotational motions, 
whereas other lines apparently do not trace the warm 
disks. For example, C'^^S traces the Keplerian motion in 
the young source IRAS 20126-1-4104 whereas it does not 
trace the central protostars at all in the sources presented 
here (Fig. [5]). In the contrary, HN^'^C shows rotational 
signatures in the HMC G29.96, but it remains completely 
undetected in the younger early-HMPO sources of our 
sample (Figs. |9]& [2]). As discussed in sections [3.2.4! and 
13.2.51 this implies that, depending on the chemical evolu- 
tion, molecules like C'^^S should be better suited for disk 
studies at very early evolutionary stages, whereas com- 
plex nitrogen-bearing molecules are promising in more 
evolved hot-core-type sources. While the chemical evo- 



aadition to tnese evolutionary etlects, o ptical aeptn is 
important for many lines. For example, ICesaroni et al.l 
(1997, 1999) have shown that CH3CN traces the rotat- 
ing structure in IRAS 20126-1-4104, whereas the same 
molec ule does not indicate any rotation in IRAS 18089- 
1732 (jBeuther et al.l l2005bf ). A hkely explanation for 
the l atter is high optical d epth of the CII3CN submm 
hnes (jBeuther et al.l[2005bl ). Again other molecules are 
excited in the accretion disks as well as the surround- 
ing envelope, causing confusion problems to disentangle 
the various physical components. In summary, getting 
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a chemical rich spectral line census like the ones pre- 
sented here shows that one can find several disk-tracing 
molecules in different sources, but it also implies that 
some previously assumed good tracers are not necessarily 
universally useful. 

The advent of broad bandpass interferometers Hke the 
SMA now fortunately allows to observe many molecular 
lines simultaneously. This way, one often finds a suit- 
able rotation-tracing molecule in an observational spec- 
tral setup. Nevertheless, one has to keep in mind the 
chemical and physical complexity in such regions, and it 
is likely that in many cases only combined modeling of 
infall, rotation and outflow will disentangle the accretion 
disk from the rest of the star-forming gas and dust core. 

4. CONCLUSION AND SUMMARY 

We compiled a sample of four massive star-forming re- 
gions in different evolutionary stages with varying lu- 
minosities that were observed in exactly the same spec- 
tral setup at high angular resolution with the SMA. We 
estimated column densities for all sources and detected 
species, and we compared the spatial distributions of the 
molecular gas. This allows us to start investigating chem- 
ical evolutionary effects also in a spatially resolved man- 
ner. Chemical modeling was conducted to explain our 
observations in more detail. 

A general result from this comparison is that many 
different physical and chemical processes are important 
to produce the complex chemical signatures wc observe. 
While some features, e.g., the non-detection of the rich 
vibrationally-torsionally excited CH3OH line forest to- 
ward the two early-HMPOs can be explained by on av- 
erage lower temperatures of the molecular gas compared 
to the more evolved HMCs, other observational charac- 
teristics require chemical evolutionary sequences caused 
by heating including grain-surface and gas phase reac- 
tions. Even other features are then better explained by 
shock-induced chemical networks. 

The rare isotopologue C^^'S is usually not detected 
right toward the main submm continuum peaks, but 
rather at the edge of the star-forming cores. This may be 
explained by temperature-selective gas-desorption pro- 
cesses and successive gas chemistry networks. Further- 
more, we find some nitrogen-bearing molecular lines to 
be only present in the HMCs, whereas they remain un- 
detected at earlier evolutionary stages. This indicates 
that the formation and excitation of many nitrogen- 
bearing molecules needs considerably higher tempera- 
tures and/or more time during the warm-up phase of the 
HMC, perhaps relating to the fact that NH3 is bonded 
within the water ice mantle. Although the statistical 
database is still too poor to set tighter constraints, these 
observations give the direction how one can use the 
presence and morphology of various molecular lines to 
identify and study different (chemical) evolutionary se- 
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quences. 

Furthermore, we discussed the observational difficulty 
to unambiguously use one or the other spectral line as 
a tracer of massive accretion disks. While some early 
spectral line candidates are discarded for such kind of 
studies by now (e.g., CN), in many other sources we 
find different lines exhibiting rotational velocity signa- 
tures. The observational feature that in most sources 
apparently only one or the other spectral line exclusively 
traces the desired structures has likely to be attributed 
to a range of effects. (1) Chemical effects, where for ex- 
ample C'^'^S may work in the youngest sources whereas 
some nitrogen-bearing molecules like HN^'^C arc better 
in typical HMCs. (2) Confusion from multiple gas com- 
ponents, mainly outflows, infall from the envelope and 
rotation. (3) High optical depth from many molecular 
lines. This implies that for future statistical studies we 
have to select spectral setups that comprise many molec- 
ular lines from various species. This way, one has good 
chances to identify for each source separately the right 
molecular tracer, and hence still draw statistically signif- 
icant conclusions. 

To advance in this field and to become more quantita- 
tive, different steps are necessary. First of all, we need to 
establish a larger database of more sources at different 
evolutionary stages, in particular even younger sources, 
as well as with varying luminosities to better characterize 
the differences and similarities. From an observational 
and technical point of view, although the presented data 
are state of the art multi-wavelength and high angular 
resolution observations, the quantitative interpretation 
is still hampered by the spatial filtering of the interfer- 
ometer. To become more quantitative, it is therefore nec- 
essary to complement such data with the missing short 
spacing information. While we have high angular resolu- 
tion in all datasets with a similar baseline coverage and 
hence similarly covered angular scales, the broad range 
of distances causes a different coverage of sampled linear 
spatial scales. Hence the missing short spacings affect 
each dataset in a different fashion which is currently the 
main limiting factor for a better quantitative interpreta- 
tion of the data. Therefore, obtaining single-dish obser- 
vations in the same spectral setup and then combining 
them with the SMA observations is a crucial step to de- 
rive more reliable column densities and from that abun- 
dances. These parameters then can be used by theorists 
to better model the chemical networks, explain the ob- 
servations and predict other suitable molecules for, e.g., 
massive disk studies. 
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Fig. 9.— Examples of rotat ion-tracing molecules : top-left: C^*S in IRAS 20126+ 4104 (ICesaroni et all [T999I . |2005|) , top-middle: 
HCOOCH3 in IRAS 18089-1732 lIBeutlier et al.ll2005bl) . top-right: H^**0 in AFGL2591 llvan der Tak et a\0m^ . bottom-left: CH3OH 
= 1 in IRAS 23151+5912 IIBeuther et al.ll2007dl) . bottom-right: HNi^C in G29.96 JBeuther et al.ll2007d '). 
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TABLE 4. Line peak intensities Speak and upper state energy levels E^/k from spectra toward peak positions of the respective 
massive star-forming regions (Fig. [2| . 
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1.6^ 
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1.6^ 
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5.9 
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" Uoubttul detection since other close vt = '2 Imcs witTi similar upper 
energy levels were not detected. 



No flux measurement possible because averaged over the given 5700 ATJ negative 
features due to missing short spacings overwhelm the positive features (see Figs. [3lfc l6t. 

Peak flux corrupted by neighboring SO2 line. 

" Only detectable with higher spatial resolution ijBeuther et al.ir2005al V 



TABLE 5 
Molecular column densities. 





Orion-KL" 


G29.96 


23151 


05358 


CH30H 


2x101^* 


4x10" 


3x10^'^ - 1x10" " 


4x10^5 - 4x10^* 


CH3CH2CN 


5x10^= 


IxlO^*^ 






CH2CHCN 


5x10^^ 


7x10^5 








2x10" 


2x101^5 


2x10" 


6xl0i3 


CH30CH3 


IxlOi** 


2x10" 






HCOOCH3 


6x10^5 


8x101*^ 






34SO 


6x10" 


IxlO^*^ 


4x101^5 




S02 


2x101^5 


3x10^'^ 


3x10^^ 




HC3N 


6x10" 


2x10^5 






HNi^c 


blend 


blend 






H2CS 


4x10" 


2x101*^ 




4x10" 


C2H5OH 


2xl0l5 


6x101*^ 






SiO 


2x10" 


4x10" 


2x10" 




CH3CN 


2xl0l5 


1x10"* 




8x101** 



" Calculated for lower average T of 200 K because of smoothing to 5700 AU resolution (Figs. |2]&;|3ll- The source size was approximated 
by half the spatial resolution. The on average lower Orion-KL column densities are likely due to the largest amount of missing flux for the 
closest source of the sample. 
From 'Beut her et_al.| l |2Q07.d ). 

From Beuther et al. (2007d) for different sub-sources. 
From Leurini et al. (2007) for different sub-sources. 

At lower temperature of 100 K, because otherwise different lines would get excited. 



